Recommendations for optimal interdisciplinary management and healthcare settings for patients with rare neurological diseases.

BACKGROUND: In 2017, the German Academy for Rare Neurological Diseases (Deutsche Akademie für Seltene Neurologische Erkrankungen; DASNE) was founded to pave the way for an optimized personalized management of patients with rare neurological diseases (RND) in all age groups. Since then a dynamic national network for rare neurological disorders has been established comprising renowned experts in neurology, pediatric neurology, (neuro-) genetics and neuroradiology. DASNE has successfully implemented case presentations and multidisciplinary discussions both at yearly symposia and monthly virtual case conferences, as well as further educational activities covering a broad spectrum of interdisciplinary expertise associated with RND. Here, we present recommendation statements for optimized personalized management of patients with RND, which have been developed and reviewed in a structured Delphi process by a group of experts. METHODS: An interdisciplinary group of 37 RND experts comprising DASNE experts, patient representatives, as well as healthcare professionals and managers was involved in the Delphi process. First, an online collection was performed of topics considered relevant for optimal patient care by the expert group. Second, a two-step Delphi process was carried out to rank the importance of the selected topics. Small interdisciplinary working groups then drafted recommendations. In two consensus meetings and one online review round these recommendations were finally consented. RESULTS: 38 statements were consented and grouped into 11 topics: health care structure, core neurological expertise and core mission, interdisciplinary team composition, diagnostics, continuous care and therapy development, case conferences, exchange / cooperation between Centers for Rare Diseases and other healthcare partners, patient advocacy group, databases, translation and health policy. CONCLUSIONS: This German interdisciplinary Delphi expert panel developed consented recommendations for optimal care of patients with RND in a structured Delphi process. These represent a basis for further developments and adjustments in the health care system to improve care for patients with RND and their families.

MRAS in coronary artery disease-Unchartered territory.

Genome-wide association studies (GWAS) have identified coronary artery disease (CAD) susceptibility locus on chromosome 3q22.3. This locus contains a cluster of several genes that includes muscle rat sarcoma virus (MRAS). Common MRAS variants are also associated with CAD causing risk factors such as hypertension, dyslipidemia, obesity, and type II diabetes. The MRAS gene is an oncogene that encodes a membrane-bound small GTPase. It is involved in a variety of signaling pathways, regulating cell differentiation and cell survival (mitogen-activated protein kinase [MAPK]/extracellular signal-regulated kinase and phosphatidylinositol 3-kinase) as well as acute phase response signaling (tumor necrosis factor [TNF] and interleukin 6 [IL6] signaling). In this review, we will summarize the role of genetic MRAS variants in the etiology of CAD and its comorbidities with the focus on tissue distribution of MRAS isoforms, cell type/tissue specificity, and mode of action of single nucleotide variants in MRAS associated complex traits. Finally, we postulate that CAD risk variants in the MRAS locus are specific to smooth muscle cells and lead to higher levels of MRAS, particularly in arterial and cardiac tissue, resulting in MAPK-dependent tissue hypertrophy or hyperplasia.

Autoantibodies against the chemokine receptor 3 predict cardiovascular risk.

BACKGROUND AND AIMS: Chronic inflammation and autoimmunity contribute to cardiovascular (CV) disease. Recently, autoantibodies (aAbs) against the CXC-motif-chemokine receptor 3 (CXCR3), a G protein-coupled receptor with a key role in atherosclerosis, have been identified. The role of anti-CXCR3 aAbs for CV risk and disease is unclear. METHODS: Anti-CXCR3 aAbs were quantified by a commercially available enzyme-linked immunosorbent assay in 5000 participants (availability: 97.1%) of the population-based Gutenberg Health Study with extensive clinical phenotyping. Regression analyses were carried out to identify determinants of anti-CXCR3 aAbs and relevance for clinical outcome (i.e. all-cause mortality, cardiac death, heart failure, and major adverse cardiac events comprising incident coronary artery disease, myocardial infarction, and cardiac death). Last, immunization with CXCR3 and passive transfer of aAbs were performed in ApoE(-/-) mice for preclinical validation. RESULTS: The analysis sample included 4195 individuals (48% female, mean age 55.5 ± 11 years) after exclusion of individuals with autoimmune disease, immunomodulatory medication, acute infection, and history of cancer. Independent of age, sex, renal function, and traditional CV risk factors, increasing concentrations of anti-CXCR3 aAbs translated into higher intima-media thickness, left ventricular mass, and N-terminal pro-B-type natriuretic peptide. Adjusted for age and sex, anti-CXCR3 aAbs above the 75th percentile predicted all-cause death [hazard ratio (HR) (95% confidence interval) 1.25 (1.02, 1.52), P = .029], driven by excess cardiac mortality [HR 2.51 (1.21, 5.22), P = .014]. A trend towards a higher risk for major adverse cardiac events [HR 1.42 (1.0, 2.0), P = .05] along with increased risk of incident heart failure [HR per standard deviation increase of anti-CXCR3 aAbs: 1.26 (1.02, 1.56), P = .03] may contribute to this observation. Targeted proteomics revealed a molecular signature of anti-CXCR3 aAbs reflecting immune cell activation and cytokine-cytokine receptor interactions associated with an ongoing T helper cell 1 response. Finally, ApoE(-/-) mice immunized against CXCR3 displayed increased anti-CXCR3 aAbs and exhibited a higher burden of atherosclerosis compared to non-immunized controls, correlating with concentrations of anti-CXCR3 aAbs in the passive transfer model. CONCLUSIONS: In individuals free of autoimmune disease, anti-CXCR3 aAbs were abundant, related to CV end-organ damage, and predicted all-cause death as well as cardiac morbidity and mortality in conjunction with the acceleration of experimental atherosclerosis.

Sex differences in the genetic and molecular mechanisms of coronary artery disease.

Sex differences in coronary artery disease (CAD) presentation, risk factors and prognosis have been widely studied. Similarly, studies on atherosclerosis have shown prominent sex differences in plaque biology. Our understanding of the underlying genetic and molecular mechanisms that drive these differences remains fragmented and largely understudied. Through reviewing genetic and epigenetic studies, we identified more than 40 sex-differential candidate genes (13 within known CAD loci) that may explain, at least in part, sex differences in vascular remodeling, lipid metabolism and endothelial dysfunction. Studies with transcriptomic and single-cell RNA sequencing data from atherosclerotic plaques highlight potential sex differences in smooth muscle cell and endothelial cell biology. Especially, phenotypic switching of smooth muscle cells seems to play a crucial role in female atherosclerosis. This matches the known sex differences in atherosclerotic phenotypes, with men being more prone to lipid-rich plaques, while women are more likely to develop fibrous plaques with endothelial dysfunction. To unravel the complex mechanisms that drive sex differences in CAD, increased statistical power and adjustments to study designs and analysis strategies are required. This entails increasing inclusion rates of women, performing well-defined sex-stratified analyses and the integration of multi-omics data.

ADAMTS-7 Modulates Atherosclerotic Plaque Formation by Degradation of TIMP-1.

BACKGROUND: The ADAMTS7 locus was genome-wide significantly associated with coronary artery disease. Lack of the ECM (extracellular matrix) protease ADAMTS-7 (A disintegrin and metalloproteinase-7) was shown to reduce atherosclerotic plaque formation. Here, we sought to identify molecular mechanisms and downstream targets of ADAMTS-7 mediating the risk of atherosclerosis. METHODS: Targets of ADAMTS-7 were identified by high-resolution mass spectrometry of atherosclerotic plaques from Apoe-/- and Apoe-/-Adamts7-/- mice. ECM proteins were identified using solubility profiling. Putative targets were validated using immunofluorescence, in vitro degradation assays, coimmunoprecipitation, and Förster resonance energy transfer-based protein-protein interaction assays. ADAMTS7 expression was measured in fibrous caps of human carotid artery plaques. RESULTS: In humans, ADAMTS7 expression was higher in caps of unstable as compared to stable carotid plaques. Compared to Apoe-/- mice, atherosclerotic aortas of Apoe-/- mice lacking Adamts-7 (Apoe-/-Adamts7-/-) contained higher protein levels of Timp-1 (tissue inhibitor of metalloprotease-1). In coimmunoprecipitation experiments, the catalytic domain of ADAMTS-7 bound to TIMP-1, which was degraded in the presence of ADAMTS-7 in vitro. ADAMTS-7 reduced the inhibitory capacity of TIMP-1 at its canonical target MMP-9 (matrix metalloprotease-9). As a downstream mechanism, we investigated collagen content in plaques of Apoe-/- and Apoe-/-Adamts7-/- mice after a Western diet. Picrosirius red staining of the aortic root revealed less collagen as a readout of higher MMP-9 activity in Apoe-/- as compared to Apoe-/- Adamts7-/- mice. To facilitate high-throughput screening for ADAMTS-7 inhibitors with the aim of decreasing TIMP-1 degradation, we designed a Förster resonance energy transfer-based assay targeting the ADAMTS-7 catalytic site. CONCLUSIONS: ADAMTS-7, which is induced in unstable atherosclerotic plaques, decreases TIMP-1 stability reducing its inhibitory effect on MMP-9, which is known to promote collagen degradation and is likewise associated with coronary artery disease. Disrupting the interaction of ADAMTS-7 and TIMP-1 might be a strategy to increase collagen content and plaque stability for the reduction of atherosclerosis-related events.

Female Gene Networks Are Expressed in Myofibroblast-Like Smooth Muscle Cells in Vulnerable Atherosclerotic Plaques.

BACKGROUND: Women presenting with coronary artery disease more often present with fibrous atherosclerotic plaques, which are currently understudied. Phenotypically modulated smooth muscle cells (SMCs) contribute to atherosclerosis in women. How these phenotypically modulated SMCs shape female versus male plaques is unknown. METHODS: Gene regulatory networks were created using RNAseq gene expression data from human carotid atherosclerotic plaques. The networks were prioritized based on sex bias, relevance for smooth muscle biology, and coronary artery disease genetic enrichment. Network expression was linked to histologically determined plaque phenotypes. In addition, their expression in plaque cell types was studied at single-cell resolution using single-cell RNAseq. Finally, their relevance for disease progression was studied in female and male Apoe-/- mice fed a Western diet for 18 and 30 weeks. RESULTS: Here, we identify multiple sex-stratified gene regulatory networks from human carotid atherosclerotic plaques. Prioritization of the female networks identified 2 main SMC gene regulatory networks in late-stage atherosclerosis. Single-cell RNA sequencing mapped these female networks to 2 SMC phenotypes: a phenotypically modulated myofibroblast-like SMC network and a contractile SMC network. The myofibroblast-like network was mostly expressed in plaques that were vulnerable in women. Finally, the mice ortholog of key driver gene MFGE8 (milk fat globule EGF and factor V/VIII domain containing) showed retained expression in advanced plaques from female mice but was downregulated in male mice during atherosclerosis progression. CONCLUSIONS: Female atherosclerosis is characterized by gene regulatory networks that are active in fibrous vulnerable plaques rich in myofibroblast-like SMCs.

The interleukin-11 receptor variant p.W307R results in craniosynostosis in humans.

Craniosynostosis is characterized by the premature fusion and ossification of one or more of the sutures of the calvaria, often resulting in abnormal features of the face and the skull. In cases in which growth of the brain supersedes available space within the skull, developmental delay or cognitive impairment can occur. A complex interplay of different cell types and multiple signaling pathways are required for correct craniofacial development. In this study, we report on two siblings with craniosynostosis and a homozygous missense pathogenic variant within the IL11RA gene (c.919 T > C; p.W307R). The patients present with craniosynostosis, exophthalmos, delayed tooth eruption, mild platybasia, and a basilar invagination. The p.W307R variant is located within the arginine-tryptophan-zipper within the D3 domain of the IL-11R, a structural element known to be important for the stability of the cytokine receptor. Expression of IL-11R-W307R in cells shows impaired maturation of the IL-11R, no transport to the cell surface and intracellular retention. Accordingly, cells stably expressing IL-11R-W307R do not respond when stimulated with IL-11, arguing for a loss-of-function mutation. In summary, the IL-11R-W307R variant, reported here for the first time to our knowledge, is most likely the causative variant underlying craniosynostosis in these patients.

pyHeart4Fish: Chamber-specific heart phenotype quantification of zebrafish in high-content screens.

Cardiovascular diseases (CVDs) are the leading cause of death. Of CVDs, congenital heart diseases are the most common congenital defects, with a prevalence of 1 in 100 live births. Despite the widespread knowledge that prenatal and postnatal drug exposure can lead to congenital abnormalities, the developmental toxicity of many FDA-approved drugs is rarely investigated. Therefore, to improve our understanding of drug side effects, we performed a high-content drug screen of 1,280 compounds using zebrafish as a model for cardiovascular analyses. Zebrafish are a well-established model for CVDs and developmental toxicity. However, flexible open-access tools to quantify cardiac phenotypes are lacking. Here, we provide pyHeart4Fish, a novel Python-based, platform-independent tool with a graphical user interface for automated quantification of cardiac chamber-specific parameters, such as heart rate (HR), contractility, arrhythmia score, and conduction score. In our study, about 10.5% of the tested drugs significantly affected HR at a concentration of 20 µM in zebrafish embryos at 2 days post-fertilization. Further, we provide insights into the effects of 13 compounds on the developing embryo, including the teratogenic effects of the steroid pregnenolone. In addition, analysis with pyHeart4Fish revealed multiple contractility defects induced by seven compounds. We also found implications for arrhythmias, such as atrioventricular block caused by chloropyramine HCl, as well as (R)-duloxetine HCl-induced atrial flutter. Taken together, our study presents a novel open-access tool for heart analysis and new data on potentially cardiotoxic compounds.

The DZHK research platform: maximisation of scientific value by enabling access to health data and biological samples collected in cardiovascular clinical studies.

The German Centre for Cardiovascular Research (DZHK) is one of the German Centres for Health Research and aims to conduct early and guideline-relevant studies to develop new therapies and diagnostics that impact the lives of people with cardiovascular disease. Therefore, DZHK members designed a collaboratively organised and integrated research platform connecting all sites and partners. The overarching objectives of the research platform are the standardisation of prospective data and biological sample collections among all studies and the development of a sustainable centrally standardised storage in compliance with general legal regulations and the FAIR principles. The main elements of the DZHK infrastructure are web-based and central units for data management, LIMS, IDMS, and transfer office, embedded in a framework consisting of the DZHK Use and Access Policy, and the Ethics and Data Protection Concept. This framework is characterised by a modular design allowing a high standardisation across all studies. For studies that require even tighter criteria additional quality levels are defined. In addition, the Public Open Data strategy is an important focus of DZHK. The DZHK operates as one legal entity holding all rights of data and biological sample usage, according to the DZHK Use and Access Policy. All DZHK studies collect a basic set of data and biosamples, accompanied by specific clinical and imaging data and biobanking. The DZHK infrastructure was constructed by scientists with the focus on the needs of scientists conducting clinical studies. Through this, the DZHK enables the interdisciplinary and multiple use of data and biological samples by scientists inside and outside the DZHK. So far, 27 DZHK studies recruited well over 11,200 participants suffering from major cardiovascular disorders such as myocardial infarction or heart failure. Currently, data and samples of five DZHK studies of the DZHK Heart Bank can be applied for.

Homozygous frameshift variant in desmoglein 2 gene causes biventricular arrhythmogenic right ventricular cardiomyopathy.

Arrhythmogenic right ventricular cardiomyopathy caused by a homozygous frameshift variant (c.2358delA) in DSG2.

Female gene networks are expressed in myofibroblast-like smooth muscle cells in vulnerable atherosclerotic plaques.

Women presenting with coronary artery disease (CAD) more often present with fibrous atherosclerotic plaques, which are currently understudied. Phenotypically modulated smooth muscle cells (SMCs) contribute to atherosclerosis in women. How these phenotypically modulated SMCs shape female versus male plaques is unknown. Here, we show sex-stratified gene regulatory networks (GRNs) from human carotid atherosclerotic tissue. Prioritization of these networks identified two main SMC GRNs in late-stage atherosclerosis. Single-cell RNA-sequencing mapped these GRNs to two SMC phenotypes: a phenotypically modulated myofibroblast-like SMC network and a contractile SMC network. The myofibroblast-like GRN was mostly expressed in plaques that were vulnerable in females. Finally, mice orthologs of the female myofibroblast-like genes showed retained expression in advanced plaques from female mice but were downregulated in male mice during atherosclerosis progression. Female atherosclerosis is driven by GRNs that promote a fibrous vulnerable plaque rich in myofibroblast-like SMCs.

Genome-wide association studies of cardiovascular disease.

Genome-wide association studies (GWAS) aim to identify common genetic variants that are associated with traits and diseases. Since 2005, more than 5,000 GWAS have been published for almost as many traits. These studies have offered insights into the loci and genes underlying phenotypic traits, have highlighted genetic correlations across traits and diseases, and are beginning to demonstrate clinical utility by identifying individuals at increased risk for common diseases. GWAS have been widely utilized across cardiovascular diseases and associated phenotypic traits, with insights facilitated by multicenter registry studies and large biobank data sets. In this review, we describe how GWAS have informed the genetic architecture of cardiovascular diseases and the insights they have provided into disease pathophysiology, using archetypal conditions for both common and rare diseases. We also describe how biobank data sets can complement disease-specific studies, particularly for rarer cardiovascular diseases, and how findings from GWAS have the potential to impact on clinical care. Finally, we discuss the outstanding challenges facing research in this field and how they can be addressed.

Autoimmune pre-disease.

Approximately 5% of the world-wide population is affected by autoimmune diseases. Overall, autoimmune diseases are still difficult to treat, impose a high burden on patients, and have a significant economic impact. Like other complex diseases, e.g., cancer, autoimmune diseases develop over several years. Decisive steps in the development of autoimmune diseases are (i) the development of autoantigen-specific lymphocytes and (often) autoantibodies and (ii) potentially clinical disease manifestation at a later stage. However, not all healthy individuals with autoantibodies develop disease manifestations. Identifying autoantibody-positive healthy individuals and monitoring and inhibiting their switch to inflammatory autoimmune disease conditions are currently in their infancy. The switch from harmless to inflammatory autoantigen-specific T and B-cell and autoantibody responses seems to be the hallmark for the decisive factor in inflammatory autoimmune disease conditions. Accordingly, biomarkers allowing us to predict this progression would have a significant impact. Several factors, such as genetics and the environment, especially diet, smoking, exposure to pollutants, infections, stress, and shift work, might influence the progression from harmless to inflammatory autoimmune conditions. To inspire research directed at defining and ultimately targeting autoimmune predisease, here, we review published evidence underlying the progression from health to autoimmune predisease and ultimately to clinically manifest inflammatory autoimmune disease, addressing the following 3 questions: (i) what is the current status, (ii) what is missing, (iii) and what are the future perspectives for defining and modulating autoimmune predisease.

Genetic and modifiable risk factors combine multiplicatively in common disease.

BACKGROUND: The joint contribution of genetic and environmental exposures to noncommunicable diseases is not well characterized. OBJECTIVES: We modeled the cumulative effects of common risk alleles and their prevalence variations with classical risk factors. METHODS: We analyzed mathematically and statistically numbers and effect sizes of established risk alleles for coronary artery disease (CAD) and other conditions. RESULTS: In UK Biobank, risk alleles counts in the lowest (175.4) and highest decile (205.7) of the distribution differed by only 16.9%, which nevertheless increased CAD prevalence 3.4-fold (p < 0.01). Irrespective of the affected gene, a single risk allele multiplied the effects of all others carried by a person, resulting in a 2.9-fold stronger effect size in the top versus the bottom decile (p < 0.01) and an exponential increase in risk (R > 0.94). Classical risk factors shifted effect sizes to the steep upslope of the logarithmic function linking risk allele numbers with CAD prevalence. Similar phenomena were observed in the Estonian Biobank and for risk alleles affecting diabetes mellitus, breast and prostate cancer. CONCLUSIONS: Alleles predisposing to common diseases can be carried safely in large numbers, but few additional ones lead to sharp risk increments. Here, we describe exponential functions by which risk alleles combine interchangeably but multiplicatively with each other and with modifiable risk factors to affect prevalence. Our data suggest that the biological systems underlying these diseases are modulated by hundreds of genes but become only fragile when a narrow window of total risk, irrespective of its genetic or environmental origins, has been passed.

Gene Sequencing Identifies Perturbation in Nitric Oxide Signaling as a Nonlipid Molecular Subtype of Coronary Artery Disease.

BACKGROUND: A key goal of precision medicine is to disaggregate common, complex diseases into discrete molecular subtypes. Rare coding variants in the low-density lipoprotein receptor gene (LDLR) are identified in 1% to 2% of coronary artery disease (CAD) patients, defining a molecular subtype with risk driven by hypercholesterolemia. METHODS: To search for additional subtypes, we compared the frequency of rare, predicted loss-of-function and damaging missense variants aggregated within a given gene in 41 081 CAD cases versus 217 115 controls. RESULTS: Rare variants in LDLR were most strongly associated with CAD, present in 1% of cases and associated with 4.4-fold increased CAD risk. A second subtype was characterized by variants in endothelial nitric oxide synthase gene (NOS3), a key enzyme regulating vascular tone, endothelial function, and platelet aggregation. A rare predicted loss-of-function or damaging missense variants in NOS3 was present in 0.6% of cases and associated with 2.42-fold increased risk of CAD (95% CI, 1.80-3.26; P=5.50×10-9). These variants were associated with higher systolic blood pressure (+3.25 mm Hg; [95% CI, 1.86-4.65]; P=5.00×10-6) and increased risk of hypertension (adjusted odds ratio 1.31; [95% CI, 1.14-1.51]; P=2.00×10-4) but not circulating cholesterol concentrations, suggesting that, beyond lipid pathways, nitric oxide synthesis is a key nonlipid driver of CAD risk. CONCLUSIONS: Beyond LDLR, we identified an additional nonlipid molecular subtype of CAD characterized by rare variants in the NOS3 gene.

Genetically Determined Reproductive Aging and Coronary Heart Disease: A Bidirectional 2-sample Mendelian Randomization.

BACKGROUND: Accelerated reproductive aging, in women indicated by early natural menopause, is associated with increased coronary heart disease (CHD) risk in observational studies. Conversely, an adverse CHD risk profile has been suggested to accelerate menopause. OBJECTIVES: To study the direction and evidence for causality of the relationship between reproductive aging and (non-)fatal CHD and CHD risk factors in a bidirectional Mendelian randomization (MR) approach, using age at natural menopause (ANM) genetic variants as a measure for genetically determined reproductive aging in women. We also studied the association of these variants with CHD risk (factors) in men. DESIGN: Two-sample MR, using both cohort data as well as summary statistics, with 4 methods: simple and weighted median-based, standard inverse-variance weighted (IVW) regression, and MR-Egger regression. PARTICIPANTS: Data from EPIC-CVD and summary statistics from UK Biobank and publicly available genome-wide association studies were pooled for the different analyses. MAIN OUTCOME MEASURES: CHD, CHD risk factors, and ANM. RESULTS: Across different methods of MR, no association was found between genetically determined reproductive aging and CHD risk in women (relative risk estimateIVW = 0.99; 95% confidence interval (CI), 0.97-1.01), or any of the CHD risk factors. Similarly, no associations were found in men. Neither did the reversed analyses show evidence for an association between CHD (risk factors) and reproductive aging. CONCLUSION: Genetically determined reproductive aging is not causally associated with CHD risk (factors) in women, nor were the genetic variants associated in men. We found no evidence for a reverse association in a combined sample of women and men.

Transcriptome-wide association study of coronary artery disease identifies novel susceptibility genes.

The majority of risk loci identified by genome-wide association studies (GWAS) are in non-coding regions, hampering their functional interpretation. Instead, transcriptome-wide association studies (TWAS) identify gene-trait associations, which can be used to prioritize candidate genes in disease-relevant tissue(s). Here, we aimed to systematically identify susceptibility genes for coronary artery disease (CAD) by TWAS. We trained prediction models of nine CAD-relevant tissues using EpiXcan based on two genetics-of-gene-expression panels, the Stockholm-Tartu Atherosclerosis Reverse Network Engineering Task (STARNET) and the Genotype-Tissue Expression (GTEx). Based on these prediction models, we imputed gene expression of respective tissues from individual-level genotype data on 37,997 CAD cases and 42,854 controls for the subsequent gene-trait association analysis. Transcriptome-wide significant association (i.e. P < 3.85e-6) was observed for 114 genes. Of these, 96 resided within previously identified GWAS risk loci and 18 were novel. Stepwise analyses were performed to study their plausibility, biological function, and pathogenicity in CAD, including analyses for colocalization, damaging mutations, pathway enrichment, phenome-wide associations with human data and expression-traits correlations using mouse data. Finally, CRISPR/Cas9-based gene knockdown of two newly identified TWAS genes, RGS19 and KPTN, in a human hepatocyte cell line resulted in reduced secretion of APOB100 and lipids in the cell culture medium. Our CAD TWAS work (i) prioritized candidate causal genes at known GWAS loci, (ii) identified 18 novel genes to be associated with CAD, and iii) suggested potential tissues and pathways of action for these TWAS CAD genes.

Cis-epistasis at the LPA locus and risk of cardiovascular diseases.

AIMS: Coronary artery disease (CAD) has a strong genetic predisposition. However, despite substantial discoveries made by genome-wide association studies (GWAS), a large proportion of heritability awaits identification. Non-additive genetic effects might be responsible for part of the unaccounted genetic variance. Here, we attempted a proof-of-concept study to identify non-additive genetic effects, namely epistatic interactions, associated with CAD. METHODS AND RESULTS: We tested for epistatic interactions in 10 CAD case-control studies and UK Biobank with focus on 8068 SNPs at 56 loci with known associations with CAD risk. We identified a SNP pair located in cis at the LPA locus, rs1800769 and rs9458001, to be jointly associated with risk for CAD [odds ratio (OR) = 1.37, P = 1.07 × 10-11], peripheral arterial disease (OR = 1.22, P = 2.32 × 10-4), aortic stenosis (OR = 1.47, P = 6.95 × 10-7), hepatic lipoprotein(a) (Lp(a)) transcript levels (beta = 0.39, P = 1.41 × 10-8), and Lp(a) serum levels (beta = 0.58, P = 8.7 × 10-32), while individual SNPs displayed no association. Further exploration of the LPA locus revealed a strong dependency of these associations on a rare variant, rs140570886, that was previously associated with Lp(a) levels. We confirmed increased CAD risk for heterozygous (relative OR = 1.46, P = 9.97 × 10-32) and individuals homozygous for the minor allele (relative OR = 1.77, P = 0.09) of rs140570886. Using forward model selection, we also show that epistatic interactions between rs140570886, rs9458001, and rs1800769 modulate the effects of the rs140570886 risk allele. CONCLUSIONS: These results demonstrate the feasibility of a large-scale knowledge-based epistasis scan and provide rare evidence of an epistatic interaction in a complex human disease. We were directed to a variant (rs140570886) influencing risk through additive genetic as well as epistatic effects. In summary, this study provides deeper insights into the genetic architecture of a locus important for cardiovascular diseases.